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Charged  particle  data  taken  by  the  S3-  3 satellite*  reported  by  Mizera  and 
Fennell  (1977)5>are  presented  as  contours~oFTfie  velocity  distribution  function 
on  a velocity- space  diagram.  This  report  focuses  on  the  analytical  technique 
used  to  interpret  the  particle  data.  Details  of  features  exhibited  by  the 
electron  and  ion  data  in  the  velocity  space  representations  are  discussed  in 
terms  of  a simple  electrostatic  acceleration  model.  The  observed  particle 
populations  are  separated  in  velocity-space  by  recognizable  demarcations 
calculated  from  the  conservation  laws  in  accordance  with  Liouville's  theorem. 
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INTRODUCTION 


Recent  charged  particle  data  acquired  by  the  S3-3  satellite,  with  good  spatial  and 
temporal  resolution,  have  been  reported  (Mizera  and  Fennell,  1977)  for  on  auroral  event  on 
August  12,  1976,  in  the  northern  auroral  zone  giving  convincing  evidence  of  the  existence 
of  a parallel  electric  field  extending  from  above  the  satellite  (7330  km)  down  to  the 
ionosphere.  The  purpose  of  this  report  is  to  elaborate  on  the  analysis  which  brings  into 
focus  the  electric  field  signature  seen  in  the  electron  and  ion  distribution  functions. 

The  features  of  the  particle  population  can  be  viewed  completely  when  a full 
satellite  spin  period  (approximately  19  seconds)  of  data  is  displayed  in  a velocity-space 
representation  as  contours  of  constant  velocity  distribution  function 

f(v)  = m2J/2E  (sec3 /km6)  , (1) 


where  J is  the  differential  particle  flux,  E is  the  energy,  m is  the  mass  and  v is  the 
velocity  of  the  ion  or  electron.  Pitch  angle  distributions  for  each  of  the  eight  energy 
channels  of  the  S3-3  electrostatic  analyzers  were  used  to  determine  the  velocity 
distribution  function  at  one-degree  intervals.  Each  velocity  distribution  function  was 
logarithmically  interpolated  to  an  arbitrary  set  of  fixed  values  of  f(v),  which  are 
logarithmically  spaced.  The  set  of  energies  corresponding  to  these  fixed  values  of  the 
distribution  function  at  each  one-degree  increment  of  pitch  angle  were  converted  to 
parallel  and  perpendicular  components  of  particle  velocity  with  respect  to  the  earth's 
magnetic  field.  This  was  done  for  both  electrons  and  ions.  Plotting  the  (v,  , vA  ) set  of 
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points  for  a constant  value  of  the  distribution  function  at  one-degree  intervals  produced  a 
contour  (Figure  1).  In  the  northern  hemisphere  the  +v^  axis  corresponds  to  particles  with 
0°  pitch  angle  precipitating  into  the  atmosphere.  The  -v(  axis  corresponds  to  particles 
with  180°  pitch  angle  coming  up  the  field  line.  Both  the  +vA  axis  and  the  -vA  axis 
correspond  to  90°  pitch  angle  particles.  The  time  sequence  of  the  acquired  data  begins  at 
and  proceeds  counterclockwise  to  360°.  Data  from  two  consecutive  half-spins  are 
displayed  on  each  plot. 

The  consequences  of  Liouville's  theorem,  applied  to  the  dynamical  trajectories  of 
charged  particles  in  phase  space  in  the  presence  of  a parallel  electric  field  has  been 
discussed  by  Whipple  (1977)  and  by  Chiu  and  Schulz  (1978).  For  reasonably  well-behaved 
electric  potentials,  expressed  as  functions  of  magnetic  field  intensity  along  a field  line, 
the  constants  of  motion  of  a charged  particle  can  be  used  to  interpret  the  data  that  have  t 
been  organized  in  the  velocity-space  diagram.  The  presence  of  a potential  drop  along  a 
magnetic  field  line  restricts  the  access  of  particle  populations  to  domains  whose 
boundaries  can  be  calculated  from  the  constants  of  motion.  For  example,  all  of  velocity- 
space  at  an  arbitrary  altitude  below  the  satellite  designated  by  the  subscript  l transforms 
into  a region  bounded  by  a hyperbola  at  satellite  altitudes.  The  equation  of  this  hyperbola 
is 

v,  2 +[1  - (B^  /B)]  vx2  - (2q/m)  |e|  (V£  -V)  = 0 , (2) 

where  B^  is  the  magnitude  of  the  earth’s  magnetic  field  at  the  altitude  denoted  by 
subscript  / (approximately  100  km  for  electrons),  B is  the  earth's  magnetic  field  at  the 
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Figure  la.  Values  of  constant  velocity- space  distribution  function  f(v)  = m J/ 
2E  (sec3-km~6)  are  plotted  for  ions  from  data  taken  on  August  12, 
1976  between  UT  12110  and  UT  12129  (Misera  and  Fennell,  1977). 
The  solid  curves  are  a family  of  hyperbolae  calculated  using 
parameters  V,  and  B(  in  equation  (2)  representing  altitudes  of 
100  km,  2000  km  and  4000  km.  The  empty  region  surrounding 
the  origin  corresponds  to  energies  below  the  threshold  of  the 
instrument.  Other  breaks  in  the  contours  are  the  result  of  data 
missing  due  to  telemetry  drop  out  or  to  particle  fluxes  which 
have  fallen  below  the  detection  threshold  of  the  instrument. 
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Figure  lb.  Values  of  constant  velocity- space  distribution  function  are 

plotted  for  electrons  corresponding  to  the  ions  shown  in  (a). 
The  solid  curves  are  the  domain  boundaries  calculated  from 
equations  (2)  and  (3)  by  using  values  of  the  local  potential  V, 
determined  by  Misera  and  Fennell  (19*77),  and  values  of  the 
parameters  Vf  and  B|  representative  of  100  km.  The  time 
sequence  of  the  data  begins  along  the  +V||  axis  corresponding 
to  particles  traveling  down  the  magnetic  field  line,  and  pro- 
ceeds through  pitch  angles  in  a counter-clockwise  direction. 
Two  consecutive  half  spin  periods  are  shown  on  the  same  plot. 
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satellite,  is  the  potential  at  the  altitude  denoted  by  l . V is  the  potential  at  the 
satellite,  q is  the  sign  of  the  charge  and  e is  the  electronic  charge.  Similarly  the  region  in 
velocity-space  which  defines  the  trajectories  of  charged  particles  that  are  turned  around 
by  the  electric  field  above  the  satellite  is  defined  by  the  interior  of  an  ellipse.  The 
equation  of  this  ellipse  is 

v,2+  [l  - (Bq/B)]  vx  2 + (2q/m)  |e|  V = 0 , (3) 

where  Bq  is  the  magnitude  of  the  earth's  magnetic  field  at  the  position  on  the  magnetic 
field  line  above  the  satellite  where  the  potential  vanishes  (arbitrarily  taken  here  to  be  at 
the  dipole  field  equator).  Figure  2,  taken  from  Chiu  and  Schulz  (1978),  shows  sche- 
matically how  the  boundaries  produced  by  equations  (2)  and  (3)  divide  velocity-space  into 
well  defined  regions.  The  regions  in  Figure  2 are  labeled  according  to  particle  origin  as 
allowed  by  access  due  to  adiabatic  motion.  For  the  purpose  of  defining  the  nomenclature 
in  Figure  2,  the  altitude  denoted  by  the  subscript/  has  been  chosen  so  that  ideally  the 
ionosphere  is  below  and  the  magnetosphere  is  above. 

The  following  discussion  will  serve  to  identify  the  various  regions  in  velocity-space 
shown  in  Figure  2.  The  region  labeled  M in  Figure  2a,  which  contains  the  origin  (v  = 0),  is 
characterized  by  a positive  ion  population  whose  mirror  points  lie  above  the  altitude 
denoted  by  / , in  the  local  hemisphere,  and  are  magnetospheric  in  nature.  In  the  -v , 
direction,  the  region  labeled  1 is  characterized  by  ions  coming  from  below  the  altitude 
denoted  by/,  that  is,  from  the  ionosphere.  The  region  along  the  +v(  axis  labeled  M,  (I)  is 
characterized  by  ions  which  are  destined  to  mirror  below  the  altitude  denoted  by  /.  These 
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ions  can  be  of  magnetospheric  origin  or  they  may  originate  in  the  ionosphere  of  the 
conjugate  hemisphere.  The  ions  in  regions  M,  (I)  have  an  energy  at  the  satellite  sufficient 
to  overcome  the  electric  potential  below.  The  ions  in  region  I have  all  been  accelerated 
by  the  electric  field  below.  The  vB  intercept  shown  in  Figure  2a  represents  the  minimum 
energy  that  ions  in  either  of  these  regions  can  have. 

The  region  above  the  hyperbola  in  Figure  2b  is  characterized  by  electrons  whose 
mirror  points  lie  above  the  altitude  denoted  by  / , whereas  the  whole  region  below  the 
hyperbola  is  characterized  by  electrons  with  mirror  points  below  thd  altitude  denoted  by 
l.  The  ellipse  is  determined  by  the  electric  potential  above  the  satellite.  Electrons  with 
velocities  interior  to  the  ellipse  do  not  have  sufficient  energy  to  overcome  the  electric 
potential  above  the  satellite  and  are  turned  around  and  accelerated  back  down  the  field 
line.  The  minimum  energy  an  electron  accelerated  through  the  electric  potential  above 
the  satellite  can  have  is  indicated  by  the  +v!t  intercept  of  the  ellipse.  The  hyperbola  and 
the  ellipse  divide  electron  velocity-space  in  Figure  2b  into  five  basic  regions.  Region  M is 
characterized  by  magnetospheric  electrons  which  mirror  above  the  ionosphere  and  which 
have  been  accelerated  by  the  electric  field  above  the  satellite.  The  region  denoted  M,  S, 
(I)  along  the  +v , axis  is  characterized  by  electrons  which  are  destined  to  mirror  in  the 
ionosphere  or  atmosphere  and  which  have  been  accelerated  above  the  satellite.  The 
region  denoted  S,  (I)  along  the  -v  u axis  is  characterized  by  electrons  whose  mirror  points 
lie  in  the  ionosphere  below.  This  electron  population  is  composed  of  high  energy 
ionospheric  electrons  and  high  energy  secondary  electrons  which  are  not  confined  by  the 
parallel  electric  field.  The  region  denoted  by  S,  I has  two  parts.  The  region  along  the 
-v  | axis  is  characterized  by  secondary  electrons  or  electrons  of  ionospheric  origin  whose 
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energy  is  not  sufficient  to  escape  the  electric  potential.  These  electrons,  seen  by  the 

satellite  as  they  travel  up  the  magnetic  field  line,  are  turned  around  between  the  satellite 

% 

and  the  equator.  They  appear  in  the  measurements  at  the  satellite  in  region  S,  I on  the 
♦v  | axis  as  they  travel  down  the  magnetic  field  line  towards  the  ionosphere.  Particular 
interest  will  be  focused  on  the  region  denoted  by  T in  Figure  2b,  which  is  bounded  by  the 
ellipse  and  the  hyperbola.  Electrons  whose  trajectories  occupy  this  region  of  velocity- 
space  are  trapped  in  the  sense  that  they  mirror  magnetically  above  the  ionosphere  below 
the  satellite  and  are  turned  around  by  the  electric  field  above  the  satellite.  Region  T is  a 
forbidden  region  in  the  sense  that  it  cannot  be  populated  by  electrons  from  the 
magnetosphere  or  from  the  ionosphere  by  invoking  a simple  one-dimensional  adiabatic 
theory.  However,  as  the  data  will  show,  a significant  electron  population  is  found  there. 
The  expression  for  the  v„  and  vA  intercepts  of  equations  (2)  and  (3)  are  appropriately 
indicated  on  Figure  2. 
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OBSERVATIONS  AND  DISCUSSION 


In  the  following  discussion  the  data  reported  by  Mizera  and  Pennell  (1977),  shown  in 
Figure  1,  will  be  examined  in  relation  to  boundaries  determined  by  equations  (2)  and  (3).  It 
will  be  ascertained  that  particle  populations  manifested  as  contours  of  velocity 
dtetribution  function  are  distinctly  different  across  these  boundaries.  To  establish  these 
boundaries,  values  of  By  , Vy  and  BQ  must  be  chosen.  Some  assumptions  have  to  be  made 
concerning  the  electric  potential  profile  along  the  magnetic  field  line  an<f  the  physical 
phenomena  which  violate  the  adiabatic  character  of  the  particle  motion.  Values  assigned 
to  these  variables  can  be  determined  empirically  or  physically  to  delineate  features  in  the 
data.  Good  delineation  of  features  in  the  electron  data  can  be  accomplished  choosii* 
values  on  a physical  basis. 

An  altitude  of  100  km  was  chosen  to  examine  the  electron  data  because  this 

altitude  is  normally  associated  with  the  calculation  of  the  electron  loss  cone.  The 

parameters  V and  Vy  were  determined  by  Mizera  and  Fennell  (1977)  by  examining  the 

energy  flux  spectra  of  electrons  and  positive  ions  at  pitch  angles  near  0°  and  110° 

respectively.  A circle  rather  than  an  ellipse  best  conforms  to  the  data  in  Figure  lb, 

implying  that  BQ/B  « 1.  The  zero  potential  point,  then,  is  a considerable  distance  above 

the  satellite.  The  equatorial  value  of  magnetic  field  strength  for  the  field  line  has  been 

chosen  for  Bq.  The  values  of  electric  potential  and  magnetic  field  strength,  consistent 

with  the  model  of  Chiu  and  Schulz  (1978),  are  V#  = 3900  volts.  V = 1000  volts.  B = M 
' * o 

gamma,  B = .058  gauss  and  By  * .579  gauss.  These  parameters  are  used  in  equations  (2) 
and  (3)  to  calculate  the  solid  curves  shown  in  Figure  1. 
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I One  is  struck  immediately  by  the  differences  in  character  of  the  particle  pop- 

ulations in  the  different  regions  of  Figure  lh.  In  particular,  the  shape  of  the  electron 
3 -6 

contour  f = 10  sec  -km  (the  contour  labeled  I)  follows  the  boundaries  corresponding  to 
region  T in  Figure  2b  very  well.  The  presence  of  a significant  population  of  trapped 
particles  suggests  an  electric  field  whose  lifetime  is  long  compared  to  the  bounce  time  of 
the  electrons.  Within  the  context  of  the  present  model,  region  T is  inaccessible  to 
charged  particles  by  adiabatic  processes.  Whipple  (1977)  has  suggested  that  population 
might  occtir  during  the  formation  of  the  parallel  electric  field. 

3 -6 

It  is  interesting  to  note  that  the  f = 10  sec  -km  electron  contour  follows  a maxi- 
mum in  the  velocity  distribution  function  at  a radius  corresponding  to  1 keV 
(v  = 1.87  x 104  km/sec)  in  the  precipitating  loss  cone  along  the  +v,  axis.  When  viewed  in 
terms  of  precipitating  flux,  this  ridge  is  seen  as  a peak  in  the  energy  flux  spectrum  at  the 
acceleration  potential,  as  predicted  by  Evans  (1974)  and  shown  in  Figure  3 of  Mizera  and 
Fennell  (1977).  Inside  the  1 keV  circle  on  the  upper  half  of  Figure  lb,  there  is  a relative 

4 

minimum  centered  at  approximately  v(  = 1.1  x 10  km/sec,  vA  = 0 in  the  region 
corresponding  to  S,  1 in  Figure  2b  along  the  +v(  axis.  The  electrons  represented  by  the 
contours  surrounding  this  minimum  are  of  local  ionospheric  origin  and  have  been  reflected 
downward  by  the  electric  field  above  the  satellite.  Contours  with  the  same  values  of 
f(v)  can  be  seen  along  the  -v  a axis  in  the  region  corresponding  to  S,  I in  Figure  2b.  These 
contours,  representing  secondary  electrons  produced  in  the  atmosphere,  merge  with  the 
loss  cone  and  are  not  closed.  The  electron  population  of  magnetospheric  origin, 
corresponding  to  region  M in  Figure  2b,  is  represented  in  Figure  lb  by  a series  of 
concentric  circles  (isotropic  except  in  the  local  loss  cone)  with  a steep  gradient.  The 
electron  populations  are  very  well  delineated  by  application  of  the  theory  of  adiabatic 
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motion,  a choice  of  potentials  inferred  from  the  data  and  selection  of  100  km  as  the  lower 
extent  of  adiabatic  motion. 

The  positive  ion  data,  determined  to  be  predominately  protons  (Mizera  and  Fennell, 
1977),  are  displayed  in  Figure  la  and  is  interpreted  by  means  of  an  analysis  similar  to  that 
used  to  interpret  the  electron  data.  The  contours  can  be  separated  into  a component 
roughly  circular  about  the  origin  (v=0)  and  a component  circular  about  a velocity 
corresponding  to  an  energy  of  1.4  keV  (v„  = -5.2  x 10  km/sec).  The  latter  represents  an 
ion  beam  coming  out  of  the  ionosphere.  Delineation  of  the  features  exhibited  in  the  ion 
contours  in  Figure  la  cannot  be  accomplished  by  a simple  choice  of  magnetic  field  and 
electric  potential  parameters  in  equation  (2).  Parameters  associated  with  two  altitudes 
which  have  a physical  significance  have  been  used  to  define  boundaries  with  which  to 
examine  the  ion  data;  altitude  100  km,  which  is  important  in  determining  the  electron  loss 
cone,  and  altitude  2000  km,  which  is  important  when  considering  proton  charge  exchange 
interactions  as  discussed  below.  Also,  the  parameters  of  a hyperbola  which  does  separate 
the  ion  beam  from  its  surroundings  is  examined  for  possible  physical  significance. 
Although  definite  conclusions  are  not  drawn  from  the  ion  data,  examination  of  the  data 
with  respect  to  the  several  curves  presented  does  provide  a perspective.  The  curve 
labeled  100  km  in  Figure  la  corresponds  to  the  atmospheric  loss  cone  as  defined  by 
equation  (2)  using  the  same  parameters  as  in  the  electron  data.  It  is  apparent  that  the 
criterion  defining  the  loss  cone  for  electrons  does  not  define  the  source  region  for  positive 
ions.  The  shape  of  the  ion  beam  contours  in  the  -v , direction  suggests  that  the  ions  are 
being  accelerated  from  source  points  located  over  a range  of  altitudes  along  the  magnetic 
field  line  above  100  km.  Chiu  and  Schulz  (1978),  using  the  CIRA  (1972)  model  of  the 
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neutral  atmosphere,  have  estimated  the  altitude  below  which  charge  exchange  becomes  an 
important  loss  mechanism  for  magnetospheric  protons  (and  above  which  protons  are 
relatively  free  of  charge  exchange  collisions  over  a full  bounce  period)  to  be  2000  km. 

It  is  interesting  to  look  at  the  ion  data  in  relation  to  the  boundary  determined  by 
positive  ions  which  mirror  at  2000  km.  The  electric  potential  representative  of  the 

2000  km  altitude  on  this  magnetic  field  line  is  determined  by  arbitrarily  assuming  a linear 

> 

potential  as  a function  of  altitude  between  the  satellite  and  100  km  altitude.  The 
hyperbola  labeled  2000  km  in  Figure  la  is  the  result  of  using  V£  = 2470  volts  and 
B ^ = 0.263  gauss  in  equation  (2). 

In  this  case  the  -v , intercept  of  the  2000  km  curve  approximates  the  apparent 
center  of  the  concentric  ion  beam  contours.  In  this  sense  the  latter  curve  is  the  best 
estimate  of  the  source  region.  That  is,  the  region  which  corresponds  to  I in  Figure  2a, 
representing  particles  of  ionospheric  origin,  lies  in  the  altitude  region  from  which  hot 
magnetospheric  precis  are  lost  due  to  charge  exchange.  The  electron  data  alone  cannot 
be  used  to  resolve  the  potential  below  the  satellite  sufficiently  to  determine  whether  the 
potential  at  100  km  is  3000  volts  or  2500  volts.  Both  values  produce  equally  satisfying 
boundaries  to  region  T.  The  inability  to  accurately  determine  the  loss  cone  boundary 
given  by  equation  (2)  is  due  to  the  fact  that  |(B^  /B)  -1 1 » 9,  requiring  a large  change 
in  (V^  -V)  to  produce  a detectable  change  in  the  vA  intercept  (ref.  Figure  2b).  The  impli- 
cation derived  from  the  ion  data  is  that  a potential  of  (approximately)  2500  volts  is 
reasonable  for  a source  altitude  of  2000  km.  Consequently  below  2000  km  neither  the 
electron  data  nor  the  collisionless  quasi-neutral  equilibrium  model  of  Chiu  and  Schulz 
(1978)  are  definitive. 
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The  ion  contours,  in  the  -v„  half  plane  of  Figure  la.  exhibit  several  points  where 
the  contours  change  direction  abruptly.  On  the  -v(  axis  there  is  a natural  location 
(-V|  = 180  km/sec)  for  the  intercept  of  the  hyperbola,  which  determines  the  value  of  the 
electric  potential  to  be  used.  A hyperbola  approximating  these  points  can  be  obtained, 
using  = 1170  volts  and  = .137  gauss  (the  4000-km  dipole  field  value),  which  might 
bf  interpreted  as  the  boundary  between  rragnctospheric  and  ionospheric  ions.  This 
hyperbola  is  shown  in  Figure  la  labeled  4000  km.  Although  this  curve  has  been  determined 
empirically,  the  electric  potential  used  is  not  an  unreasonable  physical  possibility. 

This  empirically  determined  curve  suggests  that  the  ion  source  is  extended  above 
100  km  and  may  be  significant  at  much  higher  altitudes  and  that  the  model  used  by  Chiu 
and  Schulz  (1978)  requires  refinement.  Chiu  (private  communication,  1977)  has  suggested 
that  the  particles  in  the  ion  beam  originating  in  the  ionosphere  experience  charge 
exchange  collisions  as  they  are  accelerated  up  the  magnetic  field  line.  The  low  energy  ion 
resulting  from  a charge  exchange  collision  is  accelerated  up  the  field  line  through  only  a 
fraction  of  the  potential  difference  between  the  satellite  and  the  ionosphere,  depending  on 
the  location  of  the  collision.  The  result  gives  the  appearance  of  an  extended  source. 
Considering  the  simplicity  of  the  model  used  by  Chiu  and  Schulz  (1978),  it  is  not  surprising 
that  an  exact  boundary  between  ionospheric  and  magnetospheric  ions  is  difficult  to 
determine. 

In  summary,  the  effects  of  an  electric  field  parallel  to  the  earth's  dipolar  magnetic 
field  in  the  auroral  region  has  been  oiscerned  by  organising  the  charged  particle  data  into 
contours  of  constant  velocity-epace  distribution  function.  Particle  data  taken  on  August 
12,  1976,  by  the  S3-3  satellite  have  provided  evidence  of  a parallel  electric  field  with  a 
potential  drop  of  1000  volts  above  the  satellite  and  2000  volts  below  the  satellite.  The 
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observed  particle  populations  are  separated  in  velocity-space  by  recognizable 
demarcations  calculated  from  the  conservation  laws  (which  include  the  electric  potential) 
in  accordance  with  Liouville's  theorem.  It  is  clear  that  the  high  temporal  and  spatial 
resolution  of  charged  particle  data  displayed  in  these  velocity-space  representations  is  a 
crucial  tool  for  studying  particle  dynamics  in  auroral  particle  research. 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operation*  ol  Th*  Aaroapac*  Corporation  ia  conducting 
experimental  and  theoretical  investigations  neceeeary  for  th*  evaluation  and 
applic  ation  of  ecientific  advance*  to  new  military  concept*  and  *y*t*ma.  Ver- 
•atility  and  flexibility  have  been  developed  to  a high  degree  by  th*  laboratory 
personnel  in  dealing  with  th*  many  problem*  encountered  in  th*  nation'*  rapidly 
developing  epac*  and  mietile  eyetem*.  Expertise  in  the  latest  scientific  devel- 
opment i*  vital  to  th*  accomplishment  of  tanks  related  to  these  problem*.  Th* 
laboratories  that  contribute  to  this  research  are: 

Aerophveics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lassrs. 

Chemistry  and  Physic*  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reaction*  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  th*  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics:  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  phyeics,  optical  and  acoustical 
imaging:  atmospheric  pollution;  millimeter  wav*  and  far-infrared  technology. 

Material*  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  compoeite*  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  ia 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigua-induced  fractures  in  structural  metals. 

Space  Science*  Laboratory:  Atmosphsric  and  ionospheric  physics,  radia- 
tion from  die  atmosphere,  density  and  composition  of  the  atmosphere,  aurora* 
and  airglow;  magnetospharic  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  th*  magnetosphere;  solar  physics,  studies  of  solar  mnguotic 
fields;  space  astronomy,  x-ray  astronomy;  th*  effects  of  nucloar  o^looions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  th*  effect*  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  system*. 
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